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Abstract 

It is shown that an operator can be defined in the abstract space of a random 
matrix ensemble whose matrix elements statistical distribution simulates the 
behavior of the distribution found in real physical systems. It is found that the 
key quantity that determine these distributions is commutator of the operator 
with the Hamiltonian. Application to symmetry breaking in quantum many- 
body system is discussed. 
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Ensembles of random matrix theory (RMT) (T|-|3| have had a wide application as models 
to describe statistical properties of eigenvalues and eigenf unctions of chaotic many-body 
systems. More general ensembles have also been considered [|J, in order to cover situations 
that depart from the conditions of applicability of RMT. One such class of ensembles is 
the so-called deformed Gaussian orthogonal ensemble (DGOE) |J. These "intermediate " 
ensembles are particularly useful when one wants to study the breaking of a given symmetry 
in a many-body system such as the atomic nucleus. Further, different variances of the 
deformed random matrix ensembles can interpolate any one of the three universal RMT 
ensembles, i.e., the GOE, the GUE and the GSE, and the Poisson ensemble, which would 
represent the transition between a fully chaotic situation, with no conserved quantity but 
energy, to a regular one, i.e., one with a complete set of operators that commute with 
the Hamiltonian. When discrete symmetries such as isospin, parity and time reversal are 
violated in a complex many-body environment, one relies on a description based on the 
transition between one GOE |],[7] into two coupled GOE 's in the first two cases and a GOE 
into GUE in the last case. This latter case has recently been studied by us in the case of 
disordered metal-insulator transition [||. 

The above mentioned statistical properties refer to fluctuations, around average values, 
of quantities connected to the eigenvalues and the eigenfunctions. These mean values are 
specific to the physical system being considered and, in fact, semiclassical estimates of them 
can be derived in terms of the underlying Hamiltonian ||. The fluctuations, on the other 
hand, have a quantum origin and are, in principle, universal, in the sense that the only 
information they carry about the system, is the class of underlying symmetry they belong. 
In order to compare the statistics generated by these fluctuations with the predictions of 
the ensembles, it is therefore necessary, in both cases, of eigenvalues and eigenfunctions, 
to perform some convenient rescaling of variables that eliminates their average behavior. 
In the case of the eigenvalues, they are first unfolded, which means that they are mapped 
onto new levels with a constant density made equal to one. For the eigenfunctions, the 
most convenient quantities to be statistically analyzed are not directly the components of 
the eigenfunctions, taken with respect to some particular basis, but rather matrix elements 
of a given operator. These matrix elements also show statistical fluctuations around their 
mean values |§ and they have to be subjected to some local average process that extracts 
secular variations as a function of the energy [TO]. However, what is still lacking is a direct 
comparison of matrix elements distributions with ensemble calculations. The difficulty being 
that, except in the limiting case of fully chaotic regime, nothing has been so far done in 
random matrix ensembles studies about these distributions. In the chaotic situation, with 
no conserved quantities and, as consequence, without any active selection rule, the matrix 
elements behave like components of a isotropic random vector and one expects them to 



follow the same distribution of the eigenstate components, i.e., the Porter-Thomas law |TT 
Nevertheless, this argument does not hold in intermediate situations between chaos and 
order. 

Of course, the reason why matrix elements distributions have not yet been investigated, 
in the context of matrix ensembles, is simple: it is not clear how one can define an operator 
associated to an observable in this abstract space. This is exactly the question that is being 
addressed in this letter. We want to show that it is possible to introduce some random oper- 
ator whose matrix elements would simulate, in some way, the behavior of observables found 
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in calculations and measurements performed in real physical systems. In the construction 
of this operator, we will be guided by the idea that when a system undergoes a chaos-order 
transition, a quantity that has a key role in determining the statistical behavior of the ma- 
trix elements of an operator, is the value of its commutator with the Hamiltonian.This is 
implied by the equation 

(Ei -E k ) <Ei\0\E k >=< E k | [O, H] | E k >= ihj t <E l \0\E k > (1) 

here the last equality was obtained using Schrodinger equation and assuming an operator 
with no explicit dependence on time. Eq.([l]) clearly shows that the commutator supplies the 
connection between the matrix elements of an observable and its behavior as a function of 
time. 

To see how this is reflected in the statistical distribution, suppose that we choose to look 
at the matrix elements of an observable O which becomes a conserved quantity in the regular 
regime. The distribution of these elements will undergo a transition from the Porter-Thomas 
law, at the chaotic side, to the singular distribution < Ei \ O \ E k > oc S k i, at the regular side, 
since the last term in (p]) is zero in this case. In a more general situation, as it occurs when 
reduced transition strengths are measured, what happens is that some transitions might 
become forbidden at the less chaotic regime. Then as the selection rules becomes operative, 
the first term in (|l|) vanishes for some pairs (7, k) and, as a consequence, we may say that 
we have a partial conservation of the observable which is causing the transitions. Again 
one should expect the missing transitions to cause a deviation of the statistical distributions 
from the Porter-Thomas law. 

To define the ensembles of random matrices we are going to work with, we follow the con- 
struction based on the Maximum Entropy Principle ||, that leads to a random Hamiltonian 
which can be cast in the form 



H = H + XH 1 , (2) 

where A is the parameter that controls the chaoticity of the ensemble. We will assume it 
to be defined in the domain < A < 1, in such a way that for A = 1, H — H GOE , and for 
A = 0, we have some reduced chaotic situation defined by the choice of Hq. Since we are 
specifically interested in the transitions from GOE to Poisson and from GOE to two coupled 
GOE's, the above requirements are sufficient to determine H and Hi. 



First, we consider the GOE^Poisson transition, in which case we write [[12 

N 



H — ^2 P%H Pi (3) 

i=0 

and 

N 

Hi=Y: P l H GOE P 3 (4) 

where H GOE is a N - dimensional random matrix taken from the Gaussian Orthogonal 
Ensemble, and we have introduced the projection operators Pj =| i >< % \, i = 1, . . . , N . 
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It is straightforward to verify, from the usual properties of projectors, that H = H GOE for 
A = 1 and, on the other hand, when A = 0, H becomes a diagonal matrix whose eigenvalues 
are known to the follow Poisson distribution 

Considering now the GOE— >2GOE's transition, we write |5]|| 

H = PH GOE P + QH GOE Q (5) 

and 

H x = PH GOE Q + QH GOE P (6) 

M 

where P = Pi an d Q = 1 — P. Here H is a two blocks diagonal matrix of dimensions 

i=i 

M and N — M and each block is by construction a GOE random matrix. Again, it is easily 
verified that H = H GOE for A = 1. 

Turning now to the operator, we choose it to have the form 

N 

= Y.PiH GOE P l . (7) 

i=0 

In the case of the transition towards Poisson, since O = H , we immediately derive the 
commutator relation 

[H,0] = \[H 1 ,0]. 
where the commutator on the RHS has matrix elements given by 

\TT TT 1 _ ( ttGOE ttGOE\ ttGOE (o\ 

i^o, - [tin ~ h jj ) H ij 1°) 

which obviously is a nonvanishing antisymmetric random matrix. Therefore, the statistical 
behavior of the elements of O is controlled by A and it evolves from the Porter-Thomas 
distribution to a singular delta distribution as A goes from 1 to 0. 

On the other hand, for the transition to two GOE's, it is convenient to separate the sum 
in Eq. (^) into two parts in which the first M terms define the operator Op and the other 
N — M define the operator Oq , by construction O = Op + Oq. The commutator with the 
Hamiltonian can then be written as 

[Op, PHP] + [0 Q , QHQ] + A (OpPHQ + OqQHP - QHPO P - PHQOq) (9) 

These terms have a simple interpretation. The first and the second ones are responsible for 
the transitions between states inside the blocks PHP and QHQ, respectively. On the other 
hand, the four terms inside the parenthesis cause transitions among states located inside 
different blocks. Since the latter terms in (^j are all multiplied by the parameter A, when 
A — > these transitions become forbidden. This property shows that the operator we have 
introduced is very convenient to study transition towards two coupled GOE's, a scenario 
appropriate to investigate discrete symmetry violation in complex quantum systems. 

To make our model more flexible we are going to consider, in this case of the transition 
to two GOE's, matrix elements of the generic operator 
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0' = (l-q)H + qO 



where H is given by Eq. @ and q varies between and 1. With this form O' represents 
an operator which has a conserved part. We have now a model with two parameters, the 
parameter A of the Hamiltonian which may be fixed by the fitting the eigenvalues distribution 
and the parameter q that selects the operator. For q — 1, we have just a selection rule, as q 
decreases we are introducing a localization in the matrix elements inside the selection rule. 

Following the standard procedure, we first construct, with the operator O, the normalized 
vector 

_ Q\E k > 
lGk> <E k \0*\E k > 

where | E k > with k = 1, . . . , N is an eigenvector of the Hamiltonian. From these N vectors 
we define the matrix elements 



T u =< Ei | a k > 





T k i 


2 


< 


Tki 


2 > 



which are the quantities to be statistically analyzed.. It is convenient to work with | T k i | , 
and, as mentioned above, perform a local average that extracts secular variation with the 
energies. Thus we introduce the quantities 

Vki 

where the average is done by using a Gaussian filter of variance equal to 2. It has become 
standard in the analysis of these quantities to histogram their logarithm. In the Figs. 2 and 
3, it is shown the numerical results obtained for the two transitions. 

In Fig. 1, calculations performed in the case of the GOE-Poisson transition are presented. 
The histograms of the logarithm of the matrix elements for four different values of the 
chaoticity parameter A are plotted together with three theoretical distributions: the Porter- 
Thomas, corresponding to a x 2 ty) distribution with v = 1 degree of freedom , a x 2 {v) where 
the degree of freedom v is derived directly from the "data " |T(J and finally a distribution in 



which the histograms were fitted with the superposition of two x 2 distributions. Whereas 
one x 2 distribution of the type suggested in Ref. |T0|] , is constrained to a peak always around 
zero, our calculations, however, suggest the need of a linear combination of two distribution 
as we have already proposed in Ref. M. This kind of behavior has seems to be typical of 
transitions in which the eigenstates become more and more localized. It can be understood 
as a signature of the multifractal nature of the states |13| . 



In Fig. 2, the results for the transition towards two GOE's are shown. The chaoticity 
parameter was fixed at value A = 0.032 and the parameter q, that measures the localization 
inside the selection rule, is varied. Again, the same theoretical distributions of Fig. 1 are 
also shown. At this value A, we expect to be near the case in which the two GOE's are 
completely decoupled. We see from the figure that the distributions are greatly dependent 
on the parameter q. 

We observe that in the extreme situation when we have two diagonal uncoupled blocks, 
the probability distribution can be written down explicitly. In fact, since inside each block 
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we have Porter-Thomas distributions and outside of them the matrix elements vanish, we 
have for a generic reduced strength y =\ T k[ | 2 , the distribution 



fM 1 \ 2 f^Ml ( M lV \ (M 2 \ 2 \^W 2 ( M 2 y\ M X M 2 

where Mi + M 2 = N. The case q = 0, in Fig.2, corresponds, for the above value of the 
chaoticity parameter, to an almost uncoupled situation. We see that we have only a slight 
deviation from Potter-Thomas. This means that the delta function, in the above expression, 
is practically not observed. This can explained by the fact that we are plotting the logarithm 
of the intensities. In this kind of plot the zero elements are not observed and, more than that, 
they gain a vanishing statistical weight that comes from the Jacobian of the transformation 
y — > lny. Physically this means that the selection rule is hardly detected in this kind of 
analysis if what is being considered are statistics of matrix elements of an operator without 
a conserving part. 

In conclusion, we have extended in this paper the Maximum entropy theory to the case 
of the distribution of matrix elements. Contrary to what have been suggested recently in the 
literature, we find that in the intermediate situation described by the deformed Gaussian 
ensembles, the distribution is like a sum of two x 2 ones. Our theory is quite fit to address 
the question of symmetry breaking in complex many-body systems. The application too 
isospin symmetry violation in light nuclei |14|| is underway [15 . 
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Figure Captions: 

Fig. 1 Four histograms of the logarithm of the matrix elements distributions of the 
random operator O, see text, in the case of the transition GOE — > Poisson, for the indicated 
values of the chaoticity parameter A. The calculations were done with matrices of dimension 
N = 100. 

Fig. 2 Four histograms of the logarithm of the matrix elements distributions of the 
random operator O', see text, in the case of the transition GOE — > 2GOE's, for the indicated 
values of the parameter q. The calculations were done with matrices of dimension iV = 120 
and block sizes M 1 = 40 and M 2 = 80. 
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